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A quantum information theoretic analysis of three flavor neutrino oscillations 
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Correlations exhibited by neutrino oscillations are studied via quantum information theoretic 
quantities. We show that the strongest type of entanglement, genuine multipartite entanglement, 
is persistent in the flavour changing states. We prove the existence of Bell type nonlocal features, 
in both its absolute and genuine avatars. Finally, we show that a measure of nonclassicality, dissen¬ 
sion, which is a generalization of quantum discord to the tripartite case, is nonzero for almost the 
entire range of time in the evolution of an initial electron-neutrino. Via these quantum information 
theoretic quantities capturing different aspects of quantum correlations, we elucidate the differences 
between the flavour types, shedding light on the quantum-information theoretic aspects of the weak 
force. 

PACS numbers: 


I. INTRODUCTION 


The study of correlations in quantum systems has a 
vast literature and draws its practical importance from 
potential applications to quantum technologies such as 
quantum cryptography and teleportation [T]. Recently, 
there has been a move towards extending these studies 
to systems in the domain of particle physics PH2T]. Neu¬ 
trino is a particularly interesting candidate for such a 
study. In nature, neutrinos are available in three flavors, 
viz, the electron-neutrino Vg, muon-neutrino and tau- 
neutrino Vr- Owing to their non-zero mass, they oscillate 
from one flavor to another. This has been confirmed by 
a plethora of experiments, using both natural and “man¬ 
made” neutrinos. 

Neutrino oscillations are fundamentally three flavor os¬ 
cillations. However, in some cases, it can be reduced 
to effective two flavor oscillations m- These elemen¬ 
tary particles interact only via weak interactions, conse¬ 
quently the effect of decoherence, as compared to other 
particles widely utilized for quantum information pro¬ 
cessing, is small. Numerous experiments have revealed 
interesting details of the physics of neutrinos [22ll27j . 
This paper asks what type of quantum correlations is 
persistent in the time evolution of an initial Vg or or 
Vr- It presents a systematic study of the many-faceted 
aspect of quantum correlations. Herewith, it contributes 
to the understanding how Nature processes quantum in¬ 
formation in the regime of elementary particles and, in 
particular, which aspect of quantum information is rele¬ 
vant in weak interaction processes. 

Three-flavor neutrino oscillations can be studied by 
mapping the state of the neutrino, treating it as a three- 


*Electronic address: subhashish@iitj.ac.in 
^Electronic address: akalok@iitj.ac.in 
^Electronic address: srik@poornaprajna.org 
^Electronic address: beatrix.hiesmayr@univie.ac.at 


mode system, to that of a three-qubit system [HIITT]. In 
particular, it was shown that the neutrino oscillations are 
related to the multi-mode entanglement of single-particle 
states which can be expressed in terms of flavor transi¬ 
tion probabilities. Here we take the study of such foun¬ 
dational issues further by characterizing three-flavor neu¬ 
trino oscillations by quantum correlations. This is non¬ 
trivial as quantum correlations in three-qubit systems are 
much more involved compared to their two-qubit coun¬ 
terparts. 

The present study of quantum correlations in three 
flavor neutrino oscillations can be broadly classified into 
three categories: 

• Entanglement: We study various types of the in¬ 
separability properties of the dynamics of neutrino 
oscillations via the von Neumann entropy and in 
terms of a nonlinear witness of genuine multipartite 
entanglement introduced in Ref. |28) . 


• Genuine multipartite nonlocality: Nonlocality - 
that is considered to be the strongest manifestation 
of quantum correlations - is studied in both its ab¬ 
solute and genuine tripartite facets, characterized 
by the Mermin inequalities [29] and Svetlichny in¬ 
equalities [50] . 


• Dissension: A tripartite generalization of quantum 
discord which is a measure of nonclassicality of cor¬ 
relations [5T] . 


The plan of the paper is as follows. In Sec. [H] we pro¬ 
vide a brief introduction to the phenomenology of neutri¬ 
nos and introduce the three-flavour mode entangled state 
which will be analysed using information theoretic tools. 
The core of the paper is Sec. HI where we characterize 
three flavor neutrino oscillations in terms of various facets 
of quantum correlations. We then conclude by providing 
an outlook. 
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II. THREE FLAVOUR NEUTRINO 
OSCILLATIONS 


The three flavours of neutrinos, Vf,^ and mix to 
form three mass eigenstates V 2 and 1 / 3 : 



\^T ) V ^3 / 


where U is the 3x3 PMNS (Pontecorvo-Maki-Nakagawa- 
Sakata) mixing matrix parameterized by three mixing 
angles ( 6 * 12 , O 23 and 6 * 13 ) and a CP violating phase S 
(C. ..charge conjugation, P. ..parity). Neglecting the 
CP violating phase (that has not been observed) the mix¬ 
ing matrix can be written as 

( C12C13 S12C13 S13 ^ 

—S12C23 — C12S23S13 C12C23 — S12S23S13 S23C13 

S12S23 — C12C23S13 —C12S23 — S12C23S13 C23C13 y 

( 2 ) 

where and Sy denote cos 9ij and sin 0 ^, respectively. 

Therefore, each flavor state is given by a linear super¬ 
position of the mass eigenstates, 

l^a) ^ ^ j (i^) 

k 

where a = e,/r,T; k = 1,2,3. As the massive neutrino 
states are eigenstates of the Hamiltonian with energy 
eigenvalues Ek , the time evolution of the mass eigenstates 
is given by 


Wk{t)} = e-i^>‘^\,^k) , (4) 

where \vk) are the mass eigenstates at time t = 0. 

Straightforwardly, the time evolution of flavor neutrino 
states computes to 

|^a(^)) — acKe(^) l^e) T I^m) 4“ I^t) j (b) 

with 

ac.0{t)=J2Ucke-^^^*U*p,. ( 6 ) 

k 

For example, if an electron-neutrino is produced at 
time t — 0 then its time evolution is given by 


l^e(O) aee(t) \l^e} 4“ l^/i) 4“ UeT-(t) |i^r) 5 (^) 

where 

aeeit) = |[/elPe-A^i* + |C/e2pe-S^^‘ + |[/,3|2e-t^^3*^ 

ae^(t) = + + Ue3U;3e-i^^\ 

aerit) = +Ue 2 K^e-^^^* + Ue3K^e-i^^C 


namely are in the ultrarelativistic regime, the flavor 
eigenstates are well defined in the context of quantum 
mechanics |16j . In this approximation the survival prob¬ 
abilities take the form 

= l-4^|C/„fcnc/„^f Sin2 

k>j \ / 

and the oscillation probabilities 
P V — 

J- — 

k>j \ / 

/ 

where — m^. As in the neutrino oscilla¬ 

tion experiments, the known quantity is the distance L 
between the source and the detector and not the propa¬ 
gation time t, therefore the propagation time t is replaced 
by the source and detector distance L in the above equa¬ 
tion. This is a valid approximation as all detected neu¬ 
trinos in the oscillation experiments are ultrarelativistic. 

The allowed ranges of the six oscillation parameters, 
three mixing angles and three mass squared differences, 
is obtained by a global fit to solar, atmospheric, reactor 
and accelerator neutrino data within the framework of 
three flavor neutrino oscillations. For normal ordering, 
the best fit values of three flavor oscillation parameters 
are [52] 

012 = 33.48° , 023 = 42.3° , 0i3 = 8.50° , (10) 


AtoIiC"^ 

10-5 eV2 


7.50, 


Am|i(::si Am| 2 )c^ 
10-3eV2 


2.457. (11) 


Following Ref. |16j we introduce the occupation num¬ 
ber of neutrinos by making the following correspondence: 


We) = |l)e®|0)^®|0),= |100) , 

WW ^ | 0 ),®| 1 )^®| 0 ), = | 010 ) , 

l^r) = | 0 ),® 10)^0 | 1 ),= 1001 ) . ( 12 ) 

Consequently, we can view the time evolution of a flavor 
eigenstate a = e, as a three-qubit state, i.e., 


\m) a — (^ote (t) | 100 )+a„^(t) | 010 )+a„.(t) | 001 ) . 

(13) 

Therefore, flavor oscillations can be related to time vari¬ 
ation of the tripartite entanglement of single particle 
states. 


III. STUDY OF QUANTUM INFORMATION 
THEORETIC PROPERTIES IN NEUTRINO 
OSCILLATIONS 


If we assume that the detected neutrinos have an en- Separability or the lack of separability, i.e., entangle- 
ergy of at least one MeV (the electron/positron mass), ment, is defined for a given state according to its possible 
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FIG. 1: Plot of the (normalized) flavor entropy (solid line, 
red) and the three probabilities (i^e —^ iZj,(pink, dashed) (sur¬ 
vival probability, Eq. ®), ^ (lightblue, dashed), —>■ 

(lightgreen, dashed) (oscillation probabilities, Eq. H) for 
an initial electron-neutrino state |^(t = 0))e = J^e(O) as a 
function of the distance traveled per energy L/E. 


factorization with respect to a given algebra [S3]. The 
separability problem is in general a NP-hard problem, 
and only necessary but not generally sufficient criteria 
exists to detect entanglement. For bipartite entangled 
quantum systems it suffices to ask whether the state is 
entangled or not. In the multi-partite case the problem 
is more involved since there exist different hierarchies of 
separability (defined later). We have defined the alge¬ 
bra by introducing the occupation number of the three 
flavours and our first goal is to understand the time evolu¬ 
tion of neutrino oscillation in terms of tools for classifying 
and detecting different types of entanglement. 

The next step would be to take potential measurement 
settings into account and analyse the different facets of 
correlations in the dynamic of neutrinos. In particular, 
we are interested whether there are correlations stronger 
than those predicted by any classical theory. The cor¬ 
relations are studied via two different approaches, one 
based on the dichotomy between predictions of quantum 
theory and different hidden parameter theories, and the 
other one quantifies the various information contents via 
entropies. 


A. Study of the entanglement properties 

Entanglement measures quantify how much a quan¬ 
tum state p fails to be separable. Axiomatically, it must 
be a nonnegative real function of a state which cannot 
increase under local operations and classical communi¬ 
cation (LOCC), and is zero for separable states. An 
entropic function generally quantifies the average infor¬ 
mation gain by learning about the outcome obtained by 
measuring a system. The von Neumann entropy, a quan¬ 
tum mechanical analogue of Shannon entropy, is defined 
by S(p) = —plogp and is zero for pure states and log(d) 
for the totally mixed state, where d is system dimension, 
and the log function usually refers to base 2. The entan¬ 


glement content can be computed by the entropy of the 
subsystems since the full system is pure. 

Considering the three possible partial traces of the 
three-qubit state under investigation, we obtain a con¬ 
cave function of the single-mode probabilities |ac(/ 3 (i)Pj 

i.e., with Pj := Trail but not subsystemj|lk(t))o,('k(t)|Q, 

5^flavour(|dI(t))a) — ^ ^ Pj Pj 

j=e,ii,T 

= - l«a/3WPl0gl««/3WP 

p=e,fj,,T 

- (1 - i««/3wniog(i - (14) 

p=e,fj,,T 

which we call the flavor entropy. This function is plotted 
in Fig. [^together with the survival probabilities i/e —>■ t'e, 
Eq. ([^, and the oscillation probabilities iZg —>■ r'e —>■ 

W, Eq- @- 

Since the flavor entropy <5'flavour(|^'(0)e) is nonzero for 
almost all time instances, the state is entangled. For this 
and all the following plots, we use as units, the oscil¬ 
lation period of an electron to a muon-neutrino. Since 
there are tiny changes in the behaviour in one period 
due to the existence of the third flavour, we always plot 
two periods. When the amount of all three probabilities, 
both the survival as well as the oscillation probabilities, 
become nearly equal, the flavour entropy becomes maxi¬ 
mal. Then the and Vt oscillation probabilities become 
greater than the survival probability of Ve resulting in a 
decrease in the uncertainty of the total state followed by 
an increase, when the probabilities get closer. Next, the 
uncertainty in the total state drops again and the pattern 
gets repeated. 

The entropy of all the three possible flavour entropies 
are compared in Fig. showing that for the muon- and 
tau-neutrinos the entropy is non-zero for almost all time 
instances. Compared to the electron-neutrino evolution, 
the flavor uncertainty of the other two flavours oscillates 
more rapidly and with higher amplitudes, reaching the 
maximal value more often. 

Let us now refine the picture by investigating the type 
of entanglement in neutrino oscillations. A tripartite 
pure state can be, for example written as 

|V'fe=3) = O |(('b) O 

|' 0 fe= 2 ) = \4>a) ® \4>bc), \4>b) '^\4'Ac), 

or \(j)AB) <Si \(j)c), 

\i’k=i) = Wabc, ( 15 ) 

where k gives the number of partitions dubbed the k- 
separability. If k equals the number of involved states, in 
our case fc = 3, the joint state is called fully separable, 
else it is partially separable. An important class of states 
are those that are not separable within any bi-partition, 
and are called genuinely multi-partite entangled. In gen¬ 
eral they allow for applications that outperform their 
classical counterparts, such as secret sharing [MIIS^. It 
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should be noted that since a k = 3-separable state is nec¬ 
essarily also k = 2-separable, fc-separable states have a 
nested-convex structure. 

Among the genuinely multi-partite entangled states, 
there are two subclasses known for three qubit states, 
the GHZ- and W-type of states. In Ref. [28] a general 
framework was introduced to detect and define different 
relevant multipartite entanglement subclasses and refined 
in several follow ups. In particular it has been shown to 
allow for a self-consistent classification also in a relativis¬ 
tic framework [36) . Generally, one would expect from 


a proper classification of different types of entanglement 
that for a relativistically boosted observer, which causes 
a change of the observed state, but not of the expectation 
value, that it remains in a certain entanglement class. We 
will therefore investigate this Lorentz invariant criterion, 
though let us emphasise that we do not take any rela¬ 
tivistic effects of a boosted observer into account in this 
contribution. 

The necessary criterion for a tripartite qubit state with 
one excitation (“1”) to be bipartite reads 


Qmckeip) = 2 I(001|p|010)| + 2 I(001|p|100)| + 2 |(010|p|100)| - (^(001|p|001) + (010|p|010) + (100|p|100) 

-f2v/(000|p|000) • (011|p|011) -k2v/(000|p|000) • (101|p|101) -f 2v/(000|p|000) • (110|p|110) ) <0. (16) 


If this criterion is violated the state p has no bipartite de¬ 
compositions, i.e., it is genuinely multipartite entangled. 
The positive terms are exactly the only non-zero off- 
diagonal terms of the W-state, \W) = ^{|100) -I- 1010) -I- 
|001)}, with one excitation in the computational basis, 
whereas the negative terms are only diagonal terms. Note 
that these negative terms are all zero for the IT-state in 
the given basis such that only this state obtains the max¬ 
imum value. 


Obviously, this criterion depends on the basis repre¬ 
sentation of the state p and has therefore to be opti¬ 
mized over all local unitary operations. Indeed, taking 
the “flavour basis” as the computational basis, Eq. (131, 
the unoptimized criterion becomes 


21 O^OLe WI+21 (i)l + 2| {t)\ 

— {\ciae(t)\'^ + + \o,aT(t)\'^) <0, (17) 

'-V-" 


which is not violated for all times. Gonsequently, op¬ 
timization over all local unitaries has to be taken into 
account for each time point and is plotted for an initial 
electron-, mnon- and tau-neutrino in Fig.[^ We find that 
the states at each time point are always genuine multi¬ 
partite entangled if at least two amplitudes of the state, 
Eq. (131, are non-zero, i.e., for almost all time instances. 
The results depicted in Fig.|^also prove that in the course 
of the time evolution the genuine multipartite W state 
(all amplitudes equal to ^) is reached. Hence, Nature 
exploits the maximum genuine multipartite entanglement 
in the occupation number basis. 


B. Genuine multipartite nonlocality 


We now ask the question whether in the course of the 
flavor oscillations, Bell-type nonlocality is persistent, i.e., 
there are correlations stronger than those predicted by 
any classical hidden variable theory. For that we inves¬ 
tigate the Svetlichny inequalities [30] which are a suffi¬ 
cient criteria for proving genuine tripartite nonlocality. 
In short, the idea is whether by measuring three observ¬ 
ables A, B, C and obtaining the results a, b, c, the proba¬ 
bility P{a, b, c) can be assumed to be factorizable as 

P{a,b,c) = j f(a6|A) • h(c|A) da;(A) , (18) 


where f, h are probabilities conditioned to the hidden 
variable A with the probability measure duj. The fac¬ 
torization, here chosen between the partitions A, B ver¬ 
sus C, corresponds to Bell’s locality assumption in his 
original derivation if considered for two systems. (The 
requirement of a full factorization, i.e., the additional fac¬ 
torization f(a6|A) = q(aA)-r(6A), which correspond to ab¬ 
solute locality, is explored later by the inequalities (20)). 
Then the necessary criteria for such a factorization of 
the conditioned probabilities, qubit 1 and 2 versus 3, are 
given by 


Tr {{ADC + AD'C + A'D'C - 
-^1213(7') = 

Tr {{AD'C + AD'C - A'DC' - 


A'DC')p) < 4 , 
A'D'C')p) < 4(19) 


with D = B + B' and D' = B — B'. Since we are 
not interested in a particular hidden parameter model, 
we also consider the two other bipartitions, namely 2|13 
and 3|12, and take the maximum over them. In Fig. 
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we have plotted the maximum of /“ and l’’, over all bi¬ 
partitions, for the time evolution of an initial electron-, 
muon- and tau-neutrino. In addition, each data point 
corresponds to the maximum of the optimization over 
all possible observables A,B,C. In the case of an initial 
electron-neutrino we find regions in the time evolution 
when the criterion does not detect genuine tripartite non¬ 
locality, whereas for the two other neutrino flavours we 
observe a stronger oscillating behaviour. Summing up, 
whereas genuine nonlocal correlation is largely present in 
the time evolution, there are specihc time regions when 
it vanishes. 

Requiring that for all three measurements a hidden 
parameter model should exist can be revealed by the fol¬ 
lowing set of inequalities [H] 

M“(p) = Tr {{ADC + AD'C)p) < 2 , 

M^{p) = Tr {{A'D'C - A'DC')p) < 2 , (20) 

which are connected to the Svetlichny inequality by 
/“213 = -b (see Refs. [30l |37]). These are the 
Mermin inequalities and their violation is an indicator 
of absolute nonlocality. Again we are interested in find¬ 
ing a contradiction to any hidden parameter model, thus 
we consider all bipartitions and take the maximum. The 
results are plotted in Fig. (including an optimization 
over all four arbitrary operators A, D,C, D'). For all 
times (except when the state is separable) the two in¬ 
equalities are violated when optimized over all measure¬ 
ment settings. This shows that assuming that the mode 


correlations can be simulated by an ensemble where all 
three subsystems are correlated to each other for all time 
instances is not possible. In contrast, correlations sim¬ 
ulated by a hybrid nonlocal-local ensemble, captured by 
inequalities (191, may exist for time instances close to the 
separable state, however, only for the electron-neutrino 
dynamics. 


To sum up, except for small time regions neutrino oscil¬ 
lations exhibit all the strong correlations, entanglement 
and Bell-type nonlocality, that are considered to give an 
advantage to quantum theory over classical theories for 
a number of information processing tasks. For complete¬ 
ness, in the next section we investigate the behavior of a 
measure of nonclassicality weaker than entanglement. 


C. Dissension-a measure of nonclassicality 


Classical mutual information, quantifying the informa¬ 
tion between two random variables A and B, can be 
defined by I{A:B) = H{A) — H{A\B), where H{A) = 
— log Pi is the Shannon entropy of the probabilities 
p of the outcomes of A and H{A\B) := H{A) — H{A, B) 
represents the classical conditional entropy and H {A, B) 
is the joint entropy of the pair of random variables {A, B) 
(see, e.g.. Ref. [1]). Mutual information can be general¬ 
ized for three random variables A, R, C by any of the 
following three equivalent expressions [55] 


h{A:B-.C) = H{A,B) - H{B\A) - H{A\B) - H{A\C) - H{B\C) + H{A,B\C), 

I2{A:B:C) = H{A) + H{B) + H{C) - H{A,B) - H{A,C) - H{B,C) + H{A,B,C), 

h{A-.B-.C) = H{A) + H{B) + H{C)-H{A,B)-H{A,C) + H{A\B,C) . (21) 


While the second of these expressions suggests a straight¬ 
forward quantum generalization, by replacing the Shan¬ 
non entropy by the corresponding von Neumann entropy 
S{p) = —Tr(plog/o), the first and third expressions lead 
to complications since the average conditioned entropy 
depends on the basis chosen and on the choice of the ran¬ 
dom variables A,B,C. Let us point out that in strong 
contrast to the bipartite mutual information the tripar¬ 
tite mutual information may also be negative. This is 
the case if for instance knowing the random variable C 
enhances the correlation between A and B. Following 
the concept of quantum discord [39l |40| , which quantifies 
nonclassical correlations, in Ref. |31] two measures for 
nonclassicality, called dissension, were introduced: 

Di{A:B:C) = Ji{A:B:C) - J2{A:B:C), 

D2{A-.B:C) = J3{A:B:C) - J2{A:B:C) , (22) 


where Ji are the quantum analogs of classical tripartite 
mutual information A, Eq. namely 

Ji{A:B:C) = ^(A, R) - ^(Bln^) - ^(Aln^) - S'(A|n'^) 

- S{B\U^)+S{A,B\U^), 

J2{A:B:C) = S{A) + S{B) + S{C) 

- S{A,B)- S{A,C)- S{B,C) + S{A,B,C), 
J3{A:B:C) = S{A) + S{B) + S{C) - S{A, B) - S{A,C) 

+ S'(A|n^’^) . (23) 

Here S'(A|n^) = J2kPkS{pA\u^) with 

= (nfpAsnf )/pfc and pk = Tt:{I{^pab) is the 
probability that outcome k is obtained. It is assumed 
that the basis of is chosen such as to minimize the 
uncertainty. A given state is denoted to be nonclassical 
for any departure of Di or D 2 from 0. Here Di and D 2 
deviations from zero can be associated to nonclassicality 
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FIG. 2: Plot the (not-normalized) flavour entropy S'fiavor, Eq. ( |14[ ), for the three initial flavour states (a) v^, (b) Vp,, (c)i't in 
terms of the distance traveled per energy L/E m units of the oscillation period of the two lightest neutrinos. The horizontal 
(dotted) line corresponds to the value of the IT-state. 



FIG. 3: Plot of the criterion Qoicfce detecting genuine multipartite entanglement, Eq. (161, optimized over local unities for the 
three initial flavour states (a) Ve, (b) Vp, {c)ur with respect to the distance traveled per energy L/E in units of the oscillation 
period of the two lightest neutrinos (300 data points). The criterion detects genuine multipartite entanglement if greater than 
zero and is maximal (= 1) only for the IT-state. 



FIG. 4: Plot of the Svetlichny criteria detecting genuine multipartite nonlocality, Eq. | |19[ ), optimized over possible bipartitions 
and optimized over all six different observables for the three initial flavour states (a) Ue, (b) Up, {c)vr as a function of the 
distance traveled per energy L/E in units of the oscillation period of the two lightest neutrinos (300 data points). The criterion 
detects genuine multipartite nonlocality if the value is above 4. 



FIG. 5: Plot of the maximum of and M°, Eq. (201, optimized over all involved operators for the three initial flavour states 
(a) z/e, (b) Vp, (c) Vr in terms of the distance traveled per energy L/E in units of the oscillation period of the two lightest 
neutrinos (300 data points). The criterion is above two if and only if no hidden variable model exists. 
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accessed by only one-mode or two-mode measurements, 
respectively. 

We find that J 2 is always zero, since S{A,B,C) is 
zero because the total state is pure and since <S'(^) = 
S{B,C), S{B) = S{A,C), S{C) = S{A,B), this being 
a particularity of the W class of states. Furthermore, 
any permutation of the entropy S'(7l|n^’‘^) = 0 is zero, 
since any projection onto the two-mode subspace gives a 
pure state which has zero uncertainty. Consequently, the 
relevant measures reduce in our case to 

Dii^c) = { 

S{e, ^i) - 5(e, n^) - 5(e, H") - ^(/r, U^) - H"), 

S{e, t) - S{e, n^) - ^(e, - S{t, H®) - 5(r, O'"), 

Sifi, t) - S{p, n^) - Sifi, n") - Sir, n^) - s{t, n'^)} , 
C2('I'a) = {5(e),5(M),5(r)} . (24) 

Here the three terms in the bracket of Di refer to the 
single electron-neutrino, single muon-neutrino and sin¬ 
gle tau-neutrino mode measurements, respectively. The 
three terms in the bracket of D2 refer to joint bipar¬ 
tite measurements in the muon-tau, electron-tau, muon- 
electron mode subspace, respectively (which are mini¬ 
mized to zero). 

In Fig. 1^ we plot the dissensions Di,D 2 minimized 
over all projective measurements for the time evolution 
of an initial electron-neutrino. The first notable point 
is that both measures are very sensitive to whether the 
nonclassicality is accessed by single or bipartite mea¬ 
surements and both measures are non-zero for almost 
all times. Interestingly, we find that for both measures 
min Di , min D2 and all measurement types there are time 
regions for which the value exceeds the corresponding 
value for the W-state, which has (minUi,minll 2 ) = 
(—1.738,0.918). For single measurements dissension Di 
is still considerably smaller than the values for the GHZ 
state (minZ?i = —3), in contrast to D2 where minll 2 = 
1. Moreover, a strong “twin-humped” pattern of D2 in 
the time evolution is found for joint measurements in the 
subspace of the two heavier neutrinos showing the exis¬ 
tence of the third neutrino flavour (r). 


IV. CONCLUSIONS AND OUTLOOK 

To sum up, we have computed several information the¬ 
oretic quantities detecting and classifying correlations for 
the time evolution of an initial electron-, muon- or tau- 
neutrino. We find that for almost all time instances the 
neutrino states exhibit genuine quantum features. 

We have analysed in detail the dynamics of initial neu¬ 
trino states via various types of entanglement properties, 
correlations that cannot be simulated by realistic hidden 
variable theories and non-classical correlations revealed 
by mutual information measures. In particular, dissen¬ 
sion turned out to be larger than that for the perfect 
W-state (Dicke state), for some time values, in strong 
contrast to the measures not involving measurements, 
i.e., the flavor entropy and the criterion detecting gen¬ 
uine multipartite entanglement. What physical signifi¬ 
cance this carries, if any, remains to be seen. 

Qualitatively, there are differences between an initial 
electron-neutrino and the other two neutrinos, i.e., with 
the former showing less nonclassical features when com¬ 
pared to its heavier counterparts, a point that may merit 
further scrutiny. In any case, we can conclude that foun¬ 
dational issues are more prominent in accelerator exper¬ 
iments (mainly producing muon-neutrinos) than in reac¬ 
tor experiments (mainly producing electron-neutrinos). 

The weak force, being one of the four known funda¬ 
mental forces in Nature, dominant in the flavour chang¬ 
ing process of neutrinos, reveals strong genuine quantum 
features such as also shown for weakly decaying spinless 
K-mesons [4T] or for the weakly decaying half integer 
spin hyperons [H]. The next step would be to under¬ 
stand how and whether Nature takes advantage of these 
strong quantum correlations for information processing 
in a natural setting. 
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